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Abstract 

Science Applications International Corporation is 
currently developing the Electric Propulsion 
Interactions Code, EPIC , as part of a project sponsored 
by the Space Environments and Effects Program at 
NASA Marshall Space Right Center. Now in its second 
year of development, EPIC is an interactive computer 
toolset that allows the construction of a 3-D spacecraft 
model, and the assessment of a variety of interactions 
between its subsystems and the plume from an electric 
thruster. This paper reports on the progress of EPIC 
including the recently added ability to exchange results 
the NASA Charging Analyzer Program, Nascap-2k. 
The capability greatly enhances EPICs range of 
applicability. Expansion of the toolset’s various physics 
models proceeds in parallel with the overall 
development of the software. Also presented are recent 
upgrades of the elastic scattering algorithm in the 
electric propulsion Plume Tool. These upgrades are 
motivated by the need to assess the effects of elastically 
scattered ions on the S/C for ion beam energies that 
exceed 1000 eV. Such energy levels are expected in 
future high-power (>10 kW) ion propulsion systems 
empowered by nuclear sources. 

Introduction 

More than one hundred satellites are currently operating 
in space with onboard electric propulsion, and many 
more missions, both near-Earth and deep space, are 
being considered. A critical engineering issue in the 
integration of electric propulsion (EP) systems is the 
potentially unfavorable interaction of their high-energy 
plumes with the surrounding spacecraft structure and 


diagnostic equipment. Such interactions may affect 
mission lifetime and sometimes even threaten mission 
success. As NASA considers more powerful EP 
technologies for long-duration missions the need to 
predict pertinent interactions also becomes more 
critical. Only a few programs with the goal to produce 
large-scale, 3-D, global computer tools for the 
assessment of such interactions exist - two in the 
U.S. 1,2 and one in Europe. 3 One of the two U.S. 
programs is EPIC, 

The EPIC toolset is currently developed by Science 
Applications International Corporation (SAIC) as part 
of a 2-year effort that is sponsored by the Space 
Environment and Effects (SEE) Program at the NASA 
Marshall Space Right Center (MSFC). The objective is 
to produce an integrated modeling package that can be 
used to aid the design of spacecraft with onboard 
electric propulsion. The code is built on the Module 
Integrator and Rule-based Intelligent Analytic Database 
(MIRIAD) architecture, which has formed the core of 
several NASA and Department of Defense programs. 

In general, the designer provides the following input to 
EPIC : satellite geometry and surface materials, thruster 
locations and plume parameters, case study parameters 
such as sputter yield coefficients, orbit, hours of 
thruster operation etc. The output may be in the form of 
contour plots of the plume map in space and of surface 
interactions on the 3-D spacecraft, 1-D plots along 
surfaces (e.g. erosion depth on a solar array as a 
function of distance from the thruster), integrated 
results over duration of mission (e.g. total induced 
torque in a given direction, total deposition of eroded 
material at a specific location on the S/C), list of results 
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in text format for post-processing. EPIC's various 
components have already been used to asses EP-S/C 
interactions for both the government 1,4,5 and industry. 6 

EPIC Development 

The EPIC concept, architecture and the various 
components that will comprise the code’s overall 
capability have been described in previous 
publications. 1,7 For completeness we outline here the 
main components of the toolset: 

Object Toolkit is the 3-D geometry-definition tool 
developed to facilitate interactive generation of the 
spacecraft geometry and materials. 

3-D interactions code is built based the MIRIAD 
architecture and constitutes the bulk of EPIC ' s 
plume-S/C interactions capability. 

2-D Plume Tool is used to produce the plume from 
the electric thruster. 

Ultimately, the output from of each main component of 
EPIC is displayed (graphically or otherwise) on the 
EPIC Graphics User Interface (GUI), which is also 
used to perform trade studies. The interface has four 
main panels as depicted in Fig 1. The “Parameter Input” 
and the “Auxiliary” panels also have sub panels. Each 
of these sub-panels displays parameters for a particular 
purpose. For example, the “Gridding” sub-panel allows 
for the specification of grid spacing, while the 
“Material” panel allows the designer to choose from a 
variety of default S/C materials and change the 
functional form of the sputter yield dependence on 
incident ion energy and angle. All panels are discussed 
in greater detail in reference 8. 

Both the interface and the code that integrates the 
various components of EPIC have been programmed 
using C# (C-sharp). C# is a new language developed by 
Microsoft that provides an easy-to-use object-oriented 
way of programming the windows user interface for a 
code, and that makes it extremely easy to combine 
components written in different languages using either 
the new Common Language Runtime (CLR) or the 
Component Object Model (COM) that has been the 
backbone of windows programming for years. Indeed, 
the most natural and powerful way to communicate 
with the MIRIAD data server that coordinates EPIC 
calculations is through COM. 

Object Toolkit, originally developed to create 
spacecraft surface models for the Nascap-2k plasma 
interactions code, 9 has been generalized to define 
spacecraft models for other analysis codes. In particular 


for EPIC , employment of “Special Objects” has been 
allowed to facilitate the specification of thrusters. 
“Special Objects” are objects understood by the target 
application that are included in Object Toolkit primarily 
to specify and display their location. They are therefore 
not a part of the spacecraft surface model. Each special 
object instance has an object type name ( e.g . 
‘Thruster”), an instance name (e.g. ‘Thrusterl”), and a 
set of properties that includes, at minimum, its position 
(“x”, *y\ “z”), and the color with which it is to be 
displayed (“r”, “g”, “b”). The properties can be edited 
for each special object instance in similar fashion to the 
attribute property values. 

Communication with Nascap-2k 

There are other large-scale computer tools in existence 
that have been designed to address specific areas of the 
larger “Environment Effects” problem such as 
charging, radiation, meteor impact, etc. Communication 
with such tools can greatly extend EPIC's ability to 
assess various S/C-Propulsion interactions, not 
presently within its standalone capability. For example, 
as we move towards high power S/C and/or direct drive 
systems a number of issues become increasingly 
relevant: how will high-voltage arrays alter plume 
potentials and how do electric potentials surrounding 
the S/C change in response to the plume? How are 
plume particle trajectories, especially CEX ions, 
altered? How does loss of solar array performance due 
to electron current collection (“parasitic currents”) 
increase or decrease as a result of the coupling? 

It is clearly desirable to augment EPIC ' s capabilities 
with those of other large-scale plasma and spacecraft 
simulation codes. We have demonstrated (in prototype) 
inter-application communication between EPIC and the 
Nascap-2k plasma analysis code. Nascap-2k is a fully 
three-dimensional code whose relevant capabilities 
include the calculation of space potentials consistent 
with plasma density and temperature fields as well as 
spacecraft surfaces, and tracking of charged particles in 
those potentials. Communication between EPIC and 
Nascap-2k was accomplished using the SOAP (Simple 
Object Access Protocol) protocol. 10 

In the prototype application, Nascap-2k obtains ion 
densities calculated from EPIC ' s plume model. It is 
assumed that EPIC and Nascap-2k contain the same 
spacecraft model, and agree as to its origin and 
orientation. Nascap-2k constructs a SOAP request for 
ion densities at points corresponding to the centers of 
the cells of its multiply nested cubic grid structure, and 
EPIC returns the requested data. Nascap-2k then 
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calculates potentials throughout three-dimensional 
space consistent with the specified ion densities, 
Maxwell-Boltzmann electron densities 

[n^noexp^/Tc)], and potentials on spacecraft surfaces. 

An example of electric potentials calculated by Nascap- 
2k around a S/C, in the presence of an EPIC- generated 
plume, is shown in Fig 2. The positive potentials 
generated by the high ion density are dearly visible. 
Also clearly visible is the plasma sheath near the 
spacecraft body and surrounding the thermal panel. 
Among others, these results might be used within 
Nascap-2k to (1) generate charge exchange or 
contaminant ions within the plume and track them to 
spacecraft surfaces. Such results could then be returned 
to EPICy (2) determine the influence of the spacecraft 
potentials on the plume, especially for charged 
spacecraft, and (3) study the interaction of multiple 
plumes. 

Elastic Scattering Between High-Energy Ions and 
Neutral Particles - Generalized Algorithm in EPIC 

With the revival of nuclear sources in space flight, 
many EP systems are projected to operate at high power 
levels (>10 kW) for increased specific impulse and 
thrust. Main beam ion energies in such missions can 
range in the thousands of electron-volts (>5000 eV) 
with ion flow rates exceeding a few milligrams per 
second (>3 mg/sec). In these propulsion systems the 
need to quantify the effects of high-energy ions (>300 
eV) that are produced by elastic scattering between 
main beam ions and neutral particles becomes critical. 
These ions can be scattered to angles greater than those 
associated with the divergence of the main beam and 
can therefore pose significant damage to surrounding 
S/C structures and diagnostics. 

Until now, EP/C’s Plume Tool included a semi- 
analytical formulation that had a limited range of 
applicability. 11 The work was motivated in part by the 
measurements shown in Fig 3, which were obtained in 
the plume of a 4-kW Hall thruster, the BPT4000. 12 The 
tests were conducted by the Aerospace Corporation for 
Lockheed Martin Space Systems Company. The 
BPT4000 yielded a Xenon main beam with average ion 
energy in the order of 300 eV. Here, we build on that 
work to present a more general algorithm that can be 
used to quantify the effects of these scattered ions for 
any electrostatic thruster beam. 


assumes no main beam geometry. It includes in the 
summation all nodes for which the main beam density 
and velocity are non-zero. The beam density and 
velocity profiles are those computed by the main beam 
algorithm, (b) no assumptions are made about the 
neutral density profile. It is taken to be whatever the 
Plume Tool has computed at each node and, (c) the 
machinery to compute the differential cross section at 
each angle during the integration is included, for any 
energy and any species so long the interaction 
potentials are known. Currently, only Xe+Xe + 
interaction potentials are available in EPIC. 

Determination of the Differential Cross Section. 

In Fig 3 the existence of secondary current density 
peaks in the energy spectra of the BPT4000 plume, with 
relatively high energies compared to the primary 
resonant charge exchange peak, are evident. The 
measurements were taken using a collimated retarding 
potential analyzer (RPA). For example, at an angle of 
40° the energy associated with the second maximum is 
approximately E/q=150V. These observed ion-flux 
crests show a marked energy dependence on angle. 
Because the energy dependence is roughly given by 
E b cos 2 0 where E b is the main ion beam energy and $ 

is the angle with respect to the thruster axis, these ions 
have been associated with elastic scattering (momentum 
transfer) events. To compute the flux and energy of 
these ions at any given point in space the differential 
cross section must first be determined. The previous 
algorithm in EPIC included a polynomial fit for the 
differential cross section that was applicable to 300 eV 
ions only. The formulation for its determination has 
been described elsewhere. Since it is now part of the 
generalized algorithm we outline here the basic 
approach. 

In the classical sense, the differential cross section 
associated with an atomic collision can be calculated 
from the deflection angle, 0(b,E c ,) 13 

oo 

0(b,E e ) = ;i-2b f j d - r (1) 

r VV l-(b/r) 2 -(V(r)/E e ) 

where E c is the center-of-mass collision energy, b is the 
impact parameter, r is the interatomic distance, r m is the 
classical turning point (point of nearest approach), as 
shown in the schematic below, and V(r) is the 
interatomic potential energy. 


The new algorithm differs from the previous algorithm 
in the following few areas: (a) the contribution of 
scattered ions to each node in the computational plane 
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differential cross section. The singularity corresponds 
to a minimum in the deflection function (leading to the 
so-called rainbow angle). At the high energies of 
interest in the present work (>300 eV) the deflection 
function barely exhibits a minimum as the scattering is 
almost solely governed by the repulsive part of the 
interaction potential. In such cases the classical 
approach (vs. the more rigorous quantum mechanical 
approach 14 ) is sufficiently accurate. Currently, only the 
repulsive potential function is included in EPIC. The 
coefficients that define the functional form of V(r) have 
been derived from averaged potentials by Amarouche et 
al, 13 and are also tabulated in Katz et al. n The 
derivative in equation (4) is obtained using a first order 
forward finite difference. 


The classical turning point is the largest root of the 
equation, 


r = r m 



V(r) _ 

Ec 


( 2 ) 


In the case of equal masses and a stationary target 
particle the deflection angles and energies in the 
laboratory and CM frames are simply related by, 

fr = 0/2, E c L = 2E C (6) 


which is now computed in the revised algorithm. Since 
it is impossible to distinguish between positive and 
negative deflections in the lab when measuring the 
scattering intensity at a given angle, the scattering angle 
in the center of mass frame, 0, is given by: 


where E TtL is the collision energy in laboratory frame. 
Using the relations from (6) it can be shown that the 
CM differential cross sections may be converted into 
the laboratory frame of reference using equation (7) 
below. 11 


0=0, 0 < 0 < 7t 


( 3 ) 


The elastic differential cross section is obtained from 
equation (4): 


I(0»E c ) - 


d a 
da 


CM 


(d0/db)sin0 


(4) 


where a is the collision cross section. The solid angle is 
given by d£2=2rc sin0 d0. The differential cross section 
for ions in an ion-neutral charge-exchange pair such as 
Xe + ■+■ Xe is calculated as follows: 


I i00 (6. E c ) = 4-1(0. E c ) + 4 1(71 - e. E c ) 


( 5 ) 


In EPIC , given the angle and energy for each scattering 
pair (s-point to f-node as it will be shown in the next 
section), the impact parameter is obtained by solving 
equation (1) using a Newton-Raphson method. The 
solution requires apriori knowledge of the interatomic 
potential function, V(r), which may in general have 
both attractive and repulsive contributions. The 
attractive part may lead to a singularity in the 


Ii ® 


d a _ d a 

da7 ~ d a { 


-4cosfr 


CM 


(7) 


The Xe differential cross sections for pure elastic 
scattering in the lab frame as functions of the lab angle 
fr , for E C(L =300 eV and E Ci l= 3000 eV are shown in Fig 
4. 

Determination of the density, average velocity and 
average energy of the elastically-scattered ions. We 
refer to the schematic in Fig 5 to formulate the elastic 
scattering model. It is assumed that at the scattering 
point “s’\ one scattering event occurs as the main beam 
ion flux, F b =n btS u bjS impacts one stationary neutral 
particle. Once scattered, each ion does not undergo 
another collision. Then the particle flow rate dN/dt 
(#particles N per unit time t) through the incremental 
surface dS=p 2 d£2, is given by: 


dN BT dS 
— =F h I— 
dt b p 2 


( 8 ) 


where, where dQ is the solid angle subtended by dS at 
point “s”, i = i ion (0, E c ) and p is the distance between “s” 
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and field node “f \ The particle flow rate is also equal 
to, 


— =n sf u sf dS (9) 

where, n<f and are the elastically scattered particle 

density and speed, respectively, at the field node “f\ 
Combining equations (7), (8) and (9) we obtain, 


Thus, n f (r,z) is determined by integrating the right- 

hand-side of equation (13) (using simple trapezoidal 
rule) and the volume increment from (14). 

The two components of the flux vector, F sf>r and F sf x are 
determined by, 


F -If | ( r f-yog*) 

!f ' r_|sf| | p | 


F. f ,= 


| (Zf~2 s ) 

1 |p| 


(15) 


F -n .. - 4 F bICQSft, f 

*«f n sf U sf “ 2 


( 10 ) 


Equation (10) expresses the magnitude of the flux 
vector, Ftf = |F sf | , at the field node location f(r f ,z f ,0) 

due to scattering of the main beam flux Fb by a neutral 
at the scattering point s(r s , z s , <|>) . 


It is also assumed that the distance between nodes is 
much larger than interatomic distances. Thus, we treat 
the scattering event as a single elastic collision. Since 
both momentum and kinetic energy are conserved in the 
scattering event we can write, 


and are combined with equation (13) to compute the 
average velocity components of the scattered ions as 
follows: 


Jfsf.r^o 


dV 


/Fsf. z n 0 dV 


u f,r( r >Z) = * 


u fz (r,z) = - 


(16) 


n f (r,z) ’ n f (r,z) 

As with equation (13) the integrals above are 
determined numerically. 


Finally, we compute the average energy of the scattered 
ions at each filed node using the energy flux and the 
particle flux as shown in equation (17) below, 


u rf cos 


( 11 ) 


which also yields the density of scattered ions at “f 
n jf as follows: 


F rf _ 4n b u M IcosT8 rf 4n„I 
u tf p 2 u btS cosfl rf p 2 


( 12 ) 


We now sum the contributions to “f ’ from all points “s” 
along the scattering ring defined by each scattering 
node at s(r s ,z s ,0), and by all nodes as shown in 
equation (13). 


n f (r,z)= J^n.dV 
v P 


(13) 


E f (r,z) = 


fF rf E, f n 0 dV 

V 

fcfMV 


(17) 


with E sf computed using equation (11). 


Validation. Comparisons with idealized test problems 
have been conducted for the purpose of validating the 
algorithm. The comparisons showed excellent 
agreement and are described in reference 15. Here we 
present direct comparisons with data from the BPT4000 
tests mentioned in the beginning of this section. These 
measurements were also compared with results from 
the previous numerical model in EPIC , produced by 
Katz et al n The majority of the assumptions made by 
Katz et al , have been retained in the present 
calculation, the most pertinent of which are repeated 
below. Also outlined are differences between 
assumptions made in the present and the previous 
calculations. 


where is the neutral particle density at the scattering 
node. The plume code computes the following 
incremental volume at each node: 

AV S s2nj jr s dr s dz s (14) 


Assumptions . EPIC s main beam plume model is used 
to compute the expansion of the main beam ions. 11,16 
The boundary conditions used by Katz et al at the 
thruster exit are retained in the present calculation. 
However, in contrast to the previous elastic scattering 
model, which assumed an additional 8 km/sec radial 
spread of the main-beam ion velocity at the thruster exit 
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(based on data from SPT-140 engine 17 ), no additional 
beam spread is included here. As in the previous 
calculations, the electron temperature is assumed to be 
8 eV. The computed main beam plume is shown in Fig 
6, left. The expansion of neutrals from the thruster is 
computed using an annular anode gas flow model with 
isotropic emission from the ring. 4 The hollow cathode 
neutrals are assumed to have an isotropic emission at a 
constant temperature, equal to cathode orifice 
temperature (assumed to be 538°C). The maximum 
particle density, combining thruster and cathode 
neutrals, is found to be 3.87el8 m" 3 . It is noted that in 
contrast to the elastic scattering calculations by Katz et 
a/., which assumed that all the neutral gas from the 
thruster was generated 0.1 m above the thruster 
centerline due to the simplified beam geometry 
assumed for the calculations, no such assumptions are 
required here. The chamber background density is 
assumed to be constant. The value of background 
pressure for the plume calculations is taken to be 3.3e-5 
Torr, and is higher than measured (by about a factor of 
three) but consistent with integrated current and 
performance measurements. It is noted that the 
measured values were obtained using ionization 
gauges placed at each end of the chamber (upstream 
and downstream of the thruster), close to the re-entrant 
cryopumps. 12 Using an ideal gas law the particle density 
of background neutrals is 1.06el8 m" 3 . As it will be 
shown next, also included in the comparisons are 
results for the minimum and maximum chamber 
pressures measured, 0.96e-5 and 4.56e-5 Torr, 
respectively. Ion-ion collisions have been not been 
included in the calculation. 

The computed particle density of the elastically 
scattered ions is depicted in Fig 6 (right). Also shown 
is the main beam ion density (left). In Fig 7 we compare 
results from the elastic scattering calculation for 
E f (q)fE b , with the data taken lm away from the 
thruster exit, where q is the plume half angle. It is noted 
that the values for 55° and 65° in Fig 7 are rough 
estimates from broad peaks. The average energy of the 
main beam ions, E b , is found to be 273 eV. Also 

shown for comparison in Fig 7 is the square of cos(q). 
Since the energy is averaged over the ion flux (see 
equation 17) only minor differences are exhibited at the 
different background chambers. Thus, only the case of 
3.3e-5 Torr is shown. By contrast, the computed flux of 
scattered ions lm from the thruster exit (Fig 8) is 
greatly affected by the background density, as expected. 
The result for the case of 3.3e-5 Torr (also used in the 
calculations by Katz, et at.) is in close agreement with 
the flux data. It is noted that the computed values in Fig 


8 have not accounted for the reduction of the 
elastically-scattered ion flux due to charge-exchange 
collisions, and are therefore an upper-bound estimate. 
In the calculations by Katz, et al. as much as 60% of the 
flux reaching the detector was assumed to be lost by 
charge-exchange scattering. 

Preliminary Assessment of the Effects of 
Elastic Scattering on Spacecraft Employing 
High-Power Ion Propulsion 

To obtain a first quantitative idea about possible 
concerns elastically scattered ions may pose to 
spacecraft with high power ion propulsion systems, we 
present a calculation representative of NASA’s 
Evolutionary Xenon Thruster (NEXT). NEXT is the 
next generation ion propulsion system developed by 
NASA to follow the NSTAR (NASA Solar Electric 
Propulsion Technology Applications Readiness) 
engine. 18 The system is to include (among others) a 40 
cm diameter ion thruster. Tests performed at NASA 
Glenn Research Center (GRC) using a Laboratory 
Model (LM) 40-cm engine produced performance data 
for a relatively wide power range 1. 1-6.9 kW. 18,19 We 
therefore use NEXT as a representative example to 
assess the consequences of elastically-scattered ions as 
mission requirements transition from low to high power 
ion propulsion. 

Table 1. Parameters used to compare the effects of 
elastically-scattered ions at two different power levels 
of a NEXT-class ion propulsion system. 


Parameter Case 1 Case 2 


Performance 
data (from 
Patterson, et 
al. 18 and 
Soulas, et 
al. 19 ) 

Pin (kW) 

LI 

6.9 

[Thrust eff. * Pin] 

(kW) _ _ 

0.57 

4.8 

Discharge chamber 
flowrate, m d (mg/s) 

1.92 

5.45 

Thrust, T (mN) 

49.6 

238 

Beam Current, (A) 

1.2 

3.52 

Estimated 6 / 
assumed b 
values for the 
plume 
calculations 

'Average Ion Beam 
Energy, E t (eV) 

625 

1671 

hprop. utilization, t] p 

0.855 

0.885 

“Density of neutrals at 
thruster exit, tio (m' 3 ) 

7.1el6 

16el6 

'’Beam divergence 
half-angle (deg) 

30 

30 

^Electron temp. (eV) 

1.8 

1.8 


Performance data at the minimum and maximum power 
levels, shown in Table 1 below, are used to estimate the 
required values for the main beam and elastic scattering 
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calculations. The average ion beam energy (in eV) is 
estimated in expression (18) using the thrust T and 
beam current J b from Table 1, with (q/m) being the 
charge-to-mass ratio for Xe. 


This is not unexpected since erosion due to elastic 
scattering, S^, only equations (10), (11) and (19), lead 
to the following scaling (along the panel): 


(tY 



c |E t2 cos 2 (£) + c 2 

y 

(18) 


c,E w cos 2 0&) + c 2 


*^b,2 fy n o,2 


^b.I n o,l 


( 21 ) 


The density of neutrals expelled from the thruster is 
estimated using the assumed propellant utilization, the 
discharge chamber flow rate, the thruster radius and the 
speed of neutrals. The latter is assumed to be 140 m/sec 
for both cases. The sum of the main beam and 
elastically-scattered ion density is illustrated by the 
contour line plot in Fig 9 (bottom), computed for Case 
1. To quantify the enhanced degradation of surrounding 
surfaces from beam Case 1 (625 eV beam) to Case 2 
(1671 eV beam) we place an imaginary panel covered 
by a typical spacecraft material coating (e.g. Indium Tin 
Oxide (ITO), Kapton etc.) 50 cm from the thruster also 
shown in Fig 9 (bottom). The change in erosion of the 
panel due to sputtering by the two beams, 625 eV and 
1671 eV, can be expressed by the ratio: 


•q _ Y(E } 2 ,p) I? 2 
Y(E u ,p) F tI 


09) 


where, Y(Ej,P) is the sputter yield for panel material 
(in atoms/ion) as a function of incident ion energy Ej 
(eV) and incidence angle p (rad), and F* denotes the 
incident ion flux (m ’ 2 sec* 1 ). We choose a typical sputter 
yield functional form at normal incidence that 
incorporates a linear dependence of the yield with 
energy as shown in equation (20). 

Yffj ,P = 0) = c,Ei + c 2 (atoms/ion) 

c, =0.001, c 2 =0.01 (20) 

S as a function of distance from the thruster exit along 
the panel is illustrated in Fig 9 (top). As expected, for 
z>0.5 m erosion is dominated by main beam ion flux. 
For the part of the panel that falls with the main beam, 
erosion due to main beam ions is more than five orders 
of magnitude higher than erosion due to elastically 
scattered ions. Of particular interest is the erosion 
enhancement with increased beam energy at the larger 
plume angles, where no main beam ions are reaching 
the panel. We find that the minimum value of S for 
z<0.5m is about 13, which translates to 13 times higher 
erosion for a factor of 2.7 increase in beam energy and 
8.4 times increase in the jet power (0.57 to 4.8 kW). 


For example, at z~0.5m on the panel, 0=40 deg and 
the differential cross section ratio I 2 /Ii= 0.77. 
Substitution of all remaining values from Table 1 into 
equation (21) yields 8^=13. 3 which is indeed 

approximately the computed value at this location as 
shown in Fig 9. 


Conclusions 


The development of EPIC for NASA’s SEE program 
has proceeded this year with a number of 
accomplishments towards the toolset’s final version. A 
variety of required features within each main 
component of EPIC , such as ‘Thruster” specifications 
in Object Toolkit have been successfully implemented. 
Expansion of the toolset’s range of applicability has 
also been demonstrated through a Nascap-2k 
calculation that incorporated an £77C-generated plume 
and S/C (the latter produced by Object Toolkit which is 
the geometry definition module for both codes). The 
exchange was accomplished using the SOAP (Simple 
Object Access Protocol) protocol. The need to assess 
the effects of high-energy plume ions (>300 eV) on S/C 
equipment and critical surfaces surrounding the 
propulsion system, also prompted the revision of 
EPIC's existing algorithm that previously computed 
elastically-scattered ions with energies less than 300 
eV. Results from the revised algorithm have been 
compared with lab measurements taken in the plume of 
a 4-kW Hall thruster showing good agreement with the 
angular distribution of energy and ion flux. In view of 
the need to quantify the effects of these ions, as new 
mission requirements demand EP power levels in the 
tens of kilowatts, preliminary calculations have 
compared two ion beam cases that are representative of 
the ion propulsion system NEXT. The calculations 
suggest that for plume angles beyond the maximum 
main beam divergence (i.e. >40 deg half-angle), as 
much as 13 times higher erosion of typical S/C material 
due to elastically-scattered ions may be expected as the 
power level is increased from 1.1 to 6.9 kW. The 
results reveal the need for plume measurements in high- 
power ion engines to further validate the predictive 
capability of models such as the one presented in this 
paper. 
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Fig 1. The EPIC Graphics User Interface (GUI) showing the four major panels: 3D Viewing, Parameter Input, Results 
Display, and Auxiliary. 



Fig 2. Prototype Nascap-2k calculation of potentials near spacecraft (black outlined red cells) resulting from an EPIC- 
generated plume. The plume densities were imported from EPIC using SOAP. Noted is the plasma sheath near the 
spacecraft (lower right) and near the thermal panel (center). 
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Fig 3. Collimated RPA data for the BPT-4000 at discharge power of 3kW and voltage of 300V showing the angle- 
independent, high-energy main beam peaks and the angle-dependent, elastic scattering peaks (data taken at The Aerospace 
Corporation by Jim Pollard et al . I2 ) 



Deflection Angle in Lab Frame (deg) 

300 eV 

3000 eV 


Fig 4. Classical differential cross section as a function of angle in the lab frame 1 , 4*I ion (E c ,e)*cos6, for the Xe+Xe + elastic 
scattering pair at Ec,l= 300 eV and 3000 eV. 
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Fig 5. Elastic scattering arrangement used to formulate the generalized algorithm in EPIC s plume tool. Plane defined by 
the four red comers is a portion of the plume computational plane. 



Fig 6. Computed main beam (left) and elastically-scattered ion density (right) (using the new EPIC algorithm), in the 
plume of the BPT4000 engine operating in the laboratory (P=3.3e-5 Torr). 
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Plume Angle, q (deg) 


Fig 7. Comparisons between computed energy of elastically-scattered ions E f (q)and measurements in the plume of the 
BPT400o engine operating in the lab. Plotted is the ratio of elastically-scattered ion energy over the main beam energy 
E b estimated to be approximately 273 eV. 



Plume angle, q (deg) 


Fig 8. Comparisons between computed flux of elastically-scattered ions for three chamber pressures, and measurements in 
the plume of the BPT4000 engine. 
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Fig 9. Erosion along a “test 1 ’ panel due to the plume from an ion engine operating at l.lkW and 6.9 kW. The panel is 
located 50 cm from the thruster centerline. The plume was computed using EPIC’s main beam and (recently-updated) 
elastic scattering algorithms (main beam and elastically scattered ion density contours, in m’\ for Case 1 are shown in the 
bottom). 
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